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Abstract 

We study the instability of the spectrum for a class of non-selfadjoint 
anharmonic oscillators, estimating the behavior of the instability indices (i. 
e. the norm of spectral projections) associated with the large eigenvalues of 
these oscillators. More precisely, we consider the operators 

defined on L 2 (K), with 

... . f (m + 2)7r (m + 2)7rl 

in the case where m is either equal to 1 or an even positive integer. We get 
asymptotic expansions for the instability indices, extending the results of [S] 
and [||. 
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1 Introduction 

It has been known for several years that the spectrum of a non-selfadjoint operator 
A, acting on an Hilbert space H, can be very unstable under small perturbations of 
A. In other words, unlike in the selfadjoint case, the norm of the resolvent of A near 
the spectrum can blow up much faster than the inverse distance to the spectrum. 
Equivalently, the spectrum of its perturbations A + sB, with e > and any 
B G \\B\\ < 1, is not necessarily included in the set {zeC: d(z, a (A)) < e}. 

Let A € cr(A) be an isolated eigenvalue of A, and let II \ denote the spectral 



projection associated with A. In order to understand the instability of A (in the 
above sense), we define the instability index of A as the number 



= l|n A ||, 

see [6]. 

Of course k(X) > 1 in any case, and k(X) = 1 when A is sclfadjoint. 
If Tl\ has rank 1, that is, if A is simple in the sense of the algebraic multiplicity, 
we have a convenient expression for k(A), which we shall use extensively in the 
following: if u and u* denote respectively eigenvectors of A and A* associated 
with A and A, one can easily check [4] that 

K K X ) = Ti *\T" V- 1 ) 

\{u,u*)\ 

To understand the relation between spectral instability and instability indices, we 
denote by cr £ {A) the e-pseudospectra of A, that is the family of sets, indexed by £, 

a E (A) = {z£ p(A) :\\(A- z)- 1 !! > ~\ U a(A). 



From the perturbative point of view, a E (A) can be seen as the union of the per- 
turbed spectra, in the following sense: 



a e {A)= IJ <j{A + eB). 



lieu < i 



This equivalent formulation follows from a weak version of a theorem due to Roch 
and Silbermann [2"3] . 

Instability indices are closely related to the size of e-pseudospectra around A. 
Indeed, let a E denote the connected component of <r e (A) containing A, and assume 
for simplicity that a* C\a(A) = {A} and a* is bounded. Let da* and \da* | denote 
respectively the boundary and the perimeter of a* . Then, taking into account the 
expression of the spectral projection, 



n A = -^z I {z-A)~ l dz, 

I dab 



2in 



we get [1] 



Hence, 



«(A)<~ / UA-O-'Wd^^lda*] 

Zn J dt7 x Aire 



\da*\ > 2itek{\). (1.2) 



In the finite dimensional setting at least, instability indices give a better descrip- 
tion of pseudospectra. If A € M n (C) is a diagonalizable matrix with distinct 
eigenvalues Ai, . . . , A„, Embree and Trefethen show [55] that the e-pseudospectra 
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are rather well approximated by disks of radius en(Ak) around the eigenvalues. 
More precisely, there exists e$ > such that, for all e e]0, £o[, 

|J D(\ k ,e K (X k )+0(e 2 ))ca e (A)c (J D(\ k , e K (\ k ) + 0(e 2 )). (1.3) 

In the case of an infinite dimensional space, the validity of this statement should 
be investigated, as well as the dependance on X k of the 0(e 2 ) terms. 

In the following, we will consider some anharmonic oscillators 

A(m,9) = -^ + e i9 \xr, (1.4) 

where 

|*| <mi n|— j. (1.5) 

These operators are defined on L 2 (M) by considering, first on Cq°(R), the asso- 
ciated quadratic form, which is sectorial if 9 satisfies (1.5 1, see [6]. As stated in 
[5] , its spectrum consists of a sequence of discrete simple eigenvalues, denoted in 
nondecreasing modulus order by A„ = X n (m,0), \X n \ — > +oo, and the associated 
instability indices will be denoted by K n (m, 9). 

All the spectral projections of A(m, 9) are of rank 1 (see Lemma 5 of [6]), and if 
u n denotes an eigenfunction associated with X n (m,9), then formula ( |1.1[ ) yields 

L \u n (x)\ 2 dx 

Kn (m,0) = J y n 2 )> , 1.6 
|Jr <{x)dx\ 

since in this case we have A*T = TA where T denotes the complex conjugation, 
and thus u* n = u n . 

The following pictures, which were performed by using a code inspired by [27], 
represent the boundaries of pseudospectra for some anharmonic oscillators (1.4), 
for several small values of e. The inequality (1.2) and the inclusions (1.3) show 
that the values of K n (m, 9) are related to the size, depending on e, of these small 
disks appearing around the first eigenvalues on the simulations. 
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E.-B. Davies showed in [5J and [7] that K„(m, 9) grows as n — > +00 faster than 



any power of n, that is for all m e]0, +oo[ and 9 ^ satisfying ( 1.5 1, for all a > 0, 
there exists iV = N(a, m,9) £ N such that 

Vn>N, K n (m,9)>n a . 

This statement has been improved in the case m = 2 of the harmonic oscillator 
(sometimes referred as the Davies operator), since E.-B. Davies and A. Kuijlaars 
showed [8] that k„(2, 9) grows exponentially fast as n — > +00, with an explicit rate 
o(9) : 

Um ~\ogK n {2,0) = c{6). (1.7) 

rn-+oo 

The purpose of our work is to prove that this last statement actually holds 
for the so-called complex Airy operator .4.(1,0) and for the even anharmonic os- 
cillators A(2k,9), k > 1. More precisely, we will improve the estimate by getting 
asymptotic expansions in powers of n _1 as n — > +00. 

Let us stress that such a growth of the instability indices implies that the family 
of eigenfunctions of A(m, 9) can not possess any of the "good" properties usually 
expected. It does not form a basis, neither in the Hilbert sense nor in the Riesz 
or Schauder sense (see [B], [10]). This excludes any hope of decomposing properly 
an L 2 function along the eigenspaces of the operator. 

We will prove in sectio n [4} h owever, that the eigenfunctions form complete sets of 



the L space (Theorem 1.4) 



Before stating the results of our work, let us specify some notation. Given two 
functions /, g and a real sequence (aj)j>o, we will write 



/(r) ~ 5 (t)VV-t 3 

T— f+OO * ' 

3=0 

to mean that, for all N > 1, 

f(r)=9(r) {YajT-i+OiT-"- 1 ) 




as t — ^ +00. 

We define likewise the symbol ~ . 

The following results were anounced in |18j . 
We focus on the cases m = 1 and m = 2k, k £ N* in (1.4), the first case corre- 
sponding to the complex Airy operator —j^z + e lS \x\ studied in [3] and [T3]. We 
will see in section [2] that this operator can be decomposed into the direct sum of 
its Dirichlet and Neumann realizations on M + , and we will prove the 



G 



Theorem 1.1 Let (k^(1, 9)) n >i (resp. (k„ (1, 9)) n >i) denote the instability in- 
dices of the Dirichlet (resp. Neumann) realization of the operator 

(1.8) 

on R+ ; with < |0| < 3tt/4. 

There exist real sequences (oj(0))j>i and (/3j(9))j>i such that 



k2(1,9) ~ ?«p(C(J)(n-l/4)) (l + 2«iWn-M 

v V ^ / 



(1.9) 



k£(1,<9) ~ ^M ex p(C(0)(n-3/4)) I l + ]T/3#)^ I, (1.10) 



1/4 

C(0) = 2nrnJ 2 \ sin0| and K{9) = — = 

2a/3| sin0| 

with 

J , sin 2 (20/3) o cos(g/3)sin(20/3) 

?n e = 4 / H r-g- 2 — > 0. 

V sir 9 sin 



Notice that, in the case of the complex Airy operator —j^i +ix studied in [3] and 
[14], we get 



«£(l,7r/2) ~ „, , - 



fi+E a i( 7r / 2 ) n_, 'j exp (^ (n - 1/4) y 



Regarding even anharmonic oscillators, we have a similar statement: 

Theorem 1.2 Let k E N* and 9 be such that < \9\ < (fc +^ 1)7r . If K n (2k,9) 

denotes the n-th instability index of A(2k,9) = —j^i + e l6 x 2k , then there exist 
K(2k.9) > and a real sequence (C J (2k, 6))j>i such that 



K n (2k,9) ^ ' _ V -"" 1 " i I f > /'iClk.fDu i \ , (1.11) 

as n — > +oo, with 

c M = Mgf^ > 0. (1.12) 
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where 

( tan(|fl|/(fc + l)) Y nl 

e ' fe \sm(k\Q\/{k + l))+cos{k\e\/{k + l))tM\0\l{k + l))) [ ' 

l 2(fc+l) 

<po,kix) = Im / {l-t 2k ) l ' 2 dt. (1.14) 

Jo 

Remark 1.3 In £/ie harmonic case k = 1 (Davies operator), the first term in 
) yields the Davies-Kuijlaars theorem JB$ •' 



1 / 1 \ I e ie ' A 

lim -log||II B || = c 1 (fl)=4^i , == = 2Re / 



»-»+»n V ^2cos((9/2) / I cos(6>/2) 



where f(z) = log(z + \J z 2 — 1) — z\/ z 2 — 1. 



We are also interested in the properties of the eigenfunctions of the operator 
A(m,9). The following theorem has been proved in [3] in the case of Airy op- 
erator ./I (1,(9), and in [5] in the harmonic oscillator well as for A(2k,0), 
k>2, \6\ < f . We extend the result to any operator A{2k,Q) with \Q\ < (fc +^ )x : 



Theorem 1.4 For all m = 1, 2k, k > 1, and all 9 satisfying {1.5), the eigenfunc- 
tions of A(m,9) form a complete set of the space L 2 (M). 



Theorem 1 1 . 4| and the previous estimates enable us to study the convergence of 
the operator series defining the semigroup e - tA ( m ' e ) associated with A(m, 9) when 
decomposed along the projections II n . 

The following statement extends the result of [5] in the harmonic oscillator case. 

Corollary 1.5 Let\9\ < tt/2 ande- tA{m ^ be the semigroup generated by A(m, 9), 
A n = X n (m,9) the eigenvalues ofA(m,9), andH n = H n (m,9) the associated spec- 
tral projections. 



LetT(9) = c\{9)/ cos(0/2), where c\(9) is the constant in (1.12). The series 

-t-oo 



(t) = 5>-' A »(™' e )n n (m,0) 



n=l 

is not normally convergent in cases m = 1 for any t > 0, and m — 2 for t < T{9) 
; in cases m — 2 for t > T{6), and m — 2k for any t > 0, k > 2, the series 
converges normally towards e ~ t ^( m ' e ) and, for N sufficiently large and for some 
constants C\ = C\{k,9) and C2 — C2(0), the following estimate holds: 

... s { %e Ck ^ n eM-tRe\ N ) 1 k>2 

11 -tA{m,0) 1 t _ tt \ < J VN ^ v — 

1 <NU] - ) ^exp{-2cos{9/2){t-T(9))N), k = 1, t > T 

(1-15) 

where T1<:n = III + • • •+I1jv-i denotes the projection on the first N — 1 eigenspaces. 
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The paper is organized as follows. Section [2] is devoted to the analysis of the 
complex Airy operator and the proof of Theorem |1.1| We then handle the case 
of even anharmonic oscillators in section [3| while section [4] is independent and 
dedicated to the proofs of Theorem |1.4| and Corollary |1.5[ 
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2 The complex Airy operator 

In this section we study the instability indices of the operator 

A(l,9) = -^ + e ie \x\ 

acting on L 2 (R), with \9\ < 3vr/4. The potential e i9 \ x\ is singular at x = 0, 
which will lead us to decompose the operator along its action on the even and odd 
functions of L 2 (R). 



2.1 Decomposition on the half-line 

We use the decomposition 

L 2 (R)=¥®3, (2.1) 

where *}5 and 3 denote respectively the set of even and odd functions of L 2 (IR). 
Thus we can identify a function u E L 2 (M.) to a unique couple (up,uj) G *}3 x 3 
such that u — up + uj, where 

u p( x ) = ^i u ( x ) + u (~ x )) 
is the even component of u and 

m = -(u(x) - u{-xj) 

its odd component. 

Let us now consider the restrictions of up and ui to R + , still denoted up and it/, 
and let A D (1, 0) and A (1, 9) be respectively the Dirichlet and Neumann realiza- 
tions on K + of the operator — ^p- + e l0 x. Then, a function u e L 2 (R) belongs to 
the domain of A(l,9) if and only if u P G V{A D (1,9)) and uj G V{A N {1,9)) 7 even 
and odd functions satisfying respectively the Neumann and Dirichlet conditions 



at x = 0. Hence the decomposition (2.1 1 allows to identify -4(1, 9) to the operator 

A D {1,9) 
A N (l,9 
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acting on £>(4 D (1, 9)) x V(A N (1,6)). In particular, 

a(A(l, 0)) = a(A D (l, 9)) U a(A N (l, 9)) (2.2) 

and the sequence of instability indices of .4(1, 9) is the union of those of 4 £) (1, 9) 
and 4^(1,0), labelling the eigenvalues (A n ) ra >i of 4(1,0) such that |A 3 -| < |A 3 -+i|. 
The following statement (see [3]) describes the eigenfunctions and eigenvalues of 
4.'° (1,0) and 4 (1,0). The definition and main properties of the Airy function 
will be recalled in the next paragraph. 

Proposition 2.1 

1. The spectrum of A D (1,9) is discrete and consists of the 

Af = -eT Mi , j e n, 

where fij, j > 1, are the zeroes of the Airy function Ai, /ii < and Hj+i < 

2. The spectral projection II„ associated with A® has rank 1 . 

3. For all n > 1, J R+ u^(x)dx ^ and 

|J R+ K{x)dx\ 

where 

u n (x) = Ai(n n + xe l0/3 ) (2.4) 
is an eigenfunction of A D (1,0) associated with the eigenvalue A„ ■ 

The analog result holds for 4^(1,0), replacing the zeroes fx n of the Airy function 
by its critical points v n , which satisfy v\ < and, for all n > 2, ji n < v n < fJ, n -±. 
According to (2.2 1, if (A„)„>i denote the eigenvalues of 4(1,0) (in increasing 
modulus order), then \ 2 k — e 2ie ^ 3 \^k\ and A 2 fc+i = e 2 * e / 3 | i y fc+1 1. Furthermore, 
K 2k {l,9) = k?(1,6) and K 2k+1 (l,9) = < +1 (l,0). 

Our work will be based, to a large extent, on the properties of the Airy function, 
which we recall in the next paragraph. 

2.2 The Airy function 

The Airy function Ai is the exponentially decaying solution as x — > +oo of the 
(real) Airy equation 

Ai"(x)=xAi(x), (2.5) 
which is normalised by the condition 

A*(0) ' 



ff»/3r(f) 
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Notice that Ai(x) > and Ai'(x) < for all x > 0. The Airy function can be 
extended in the whole complex plane as an holomorphic function, and its behavior 
at infinity is described by asymptotic expansions (see [1]), as well as the behavior 
of its derivative, zeroes and critical points. 

In the following, we choose the principal value of the square root. The asymptotic 
properties of the Airy function are summarized in the following proposition, and 
will be used in the next paragraph to compute the instability indices k^(1,6). 

Proposition 2.2 There exist Ck, dk, ak, bk, k > 1, such that 
1. For any a g]0, tt], 



for | arg z\ < ir — a ; 



fe=i 



2. 

Ai{~z) 



^ + \) (i + D- 1 ) fc ^(|* 8/a r afc ) 



COs(|z 3 / 2 + I) g(-l)fe C2fe+1 (^3/2 r2fe -l ) {2 7) 



3 ' ^-r-v 3 

fe=0 



/or | arg z\ < \ ; 
3. For any a €]0,7r], 



1/4 / +OQ r, \ 

Ai'{z) ~ -* e -8^ l + ^C-l)**^)-* (2-8) 



/or | arg z\ < ir — a 

4- 



Aif(-z) - (cos(|z 3 / 2 + f) ( l + f (-lfrf 2fe( | 2 3/ 2) 



-2fc 



sin (I- 3/2 + l)E(- 1 ) fe ^ + i(^ 3/2 )- 2 ^ 1 ') (2-9) 
fc=0 '* 



/or | arg2r| < ^ ; 

5. 



3?r \ 2/3 / ±50 



Kl ~ -^(n-1/4) l + ^(-l) fc+1 a fc (n-l/4)- 2 M ; (2.10) 



fc=i 



11 



6. 



n— > + og 



3?r 



\ 2/3 / +00 \ 

(n-3/4)J h+^(-l) fe fe fc (n-3/4)- 2 M ; 



(2.11) 




Figure 3: Graph of the Airy function Ai on 

We also recall the following useful lemma: 
Lemma 2.3 The function 

F : x i-> xAi 2 (x) - Ai' 2 (x) 

is a primitive for Ai 2 . 
Proof: Use equation (]2.5h : 



F'(:r) = Ai 2 (a;) + 2a;Ai(a;)v4i'(x)-2 J 4i"(2;)Ai'(a;) 
= Ai 2 (a;) + 2Ai'(x)(xAi(x) - At" (a;)). 



(2.12) 



We are now ready to compute the asymptotics of the instability indices for 
A D (l,6) and A N {1,6). 
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2.3 Instability indices 

We restrict our study to the Dirichlet realization ^^(l,^) of the complex Airy 
operator, the case of the Neumann realization being similar, if we take care of 
replacing the zeroes \i n by the critical points v n , and the expansion (2.10) by 
(2.11 ). We also assume 9 > without loss of generality (in the case 9 < 0, all the 
computations and statements below hold with \9\ instead of 9). 

Let u n (x) = Ai(fj, n + e^^x) be the n-th eigenfunction of A D (1,9), and let 



u n (t) 2 dt 



In = / \u n {t)\ dt, J n 



be respectively the numerator and denominator of expression (2.3). We will de- 
termine asymptotic expansions as n — > +oo for both integrals. 

Step 1: the /„ term 

Let e > be fixed and small enough (see further specifications on the size of e). 
Then the decomposition 

Rn,e, (2.13) 



In In.e 



where 



In.e • 



\un(t)\ 2 dt and R n . E 



e\ n n\ 



\u n (t)\ 2 dt, 



will allow us to apply (2.6) and (2.7) uniformly in each integral R n ,s and I n . e 
respectively. 



Step l.a. We first handle the I n , E term. According to (2.6), there exists a 
sequence (c' k )k>i such that, for all N > 1, 



In.i 



1 



s|/*n| 



47r|/i n + te <e / 3 |V2 exp ^ 3 



-Re 



( M „+te^ 3 ) 3 / 2 



N 



1 + ]T c' fe (Mn + ie^/ 3 )- 3fc / 2 + r JV ( M „ + te lS / 3 ) 



k=l 



dt 



where 



MOI = 0(\Z\- 3N/2 ), |£|->+oo, 
since \fx n + te l8 ^ 3 \ — > +oo as n — > +oo uniformly with respect to t & [e\fj, n \, +oo[ : 

Vf > 0, \fj, n + te l9/3 \ > |tsin0/3| > sin6>/3| -> +oo, n -> +oo. 

Putting 



we then get 



h = h n :=\n n \ 3/2 and t = h 2/3 s = \fi n \s, 



(2.14) 
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where 

1 r + °° 

ie{k) = 4nh^ a e {s,h)e-^d S , (2.15) 



with 

a e (s,h) hZo - 

and 



1 / + °° 

\ A.: — 1 



= ~B^-l + *e W/3 ) 3/a , 
the functions £fc(s) satisfying, for all A; > 1, 

|4( S )|=0((l + S )- 3fe / 2 ). 



Now the right-hand side of (2.151 is in the right form to apply the Laplace method 
as h — > 0. Let us check that the function ipg has a unique non-degenerate critical 
point sg in the interval [e,+oo[, and that ip'g(sg) > 0. We have 

^00 - 2Re ( e «/3(_i + se «/3)i/2) 

= 2| - 1 + se^/ 3 ! 1 ^ C o S ^ + i arg (-l + se ie / 3 )) . 

Since arg(— 1 + se' 61 / 3 ) €]9/3, tt[ and 6* < 37r/4, ^ vanishes if and only if 



which yields 



S 1 , ;ii;o. 7T 

- + -arg(-l + S e l9/3 ) = -, 



26> s e sin(0/3) 
tan — = 



3 1 - sg cos(0/3) ' 

ll6IlC6 

tan(2fl/3) = sin(2fl/3) 

' 9 sin(0/3) + cos((9/3) tan(26»/3) sin0 ' [ ' 1 

and wee see that sg > for < 6 < 37r/4. 
On the other hand, one can easily check that 

ip' e \se) = m~ 1/2 sm9 > (2.17) 

where rag := \ — 1 + sge l6, / 3 | is the constant appearing in Theorem [LTj Finally, 
the principal term of a(sg, /i) is equal to 

a{s g ,0) = rrig 1/2 ^ 0, 

and 

I \ 4 3 / 2 • a 
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We choose e < sg and M > sg fixed, so that fg(e) > tpg(sg) and (pg(M) > 
ipg(sg). We write 

i £ {h)=i eM {h) + R M (h), (2.18) 

with 

{ IeMh) = I ^J^a e (s,h)e-^ds 
\ Ru(h) = ^J^ae(s,h)e-i^ds. 

The Laplace method applies to I £ ,M{h) and yields, according to the computations 
above, 



Ie,M( h ) 



/l-S-0 



2 y/2jrm V 4 \/sin0 



exp(— m^sinfl) [ 1 + ^ A k (6)h k ) (2.19) 



fe=i 



for a real sequence (Ak(6))k>i- 

Regarding R,M(h), since <pg is non-decreasing on [M, +oo[, we have 



< K'h- X l z e-^ M \ K'>0, 

as soon as /i < 1, and of course <pg(M) > ipg(sg). 
Using pTl8l ), |2l9| and ( jOo] ), we get 



+°° 2 v / 2^my 4 Vsin^ V 3/l 



+oo 



fc=l 



(2.20) 



exp — m* /2 sin0 l + ^A fe (#)/i fe . (2.21) 



Remark 2.4 In the case of the operator —j^i + ix, which is studied in JSj/ and 
the critical point of <Pn/2 * s reached when 



arg(— 1 + se 



i0/3\ 



2n 



We then get the following values: s^/ 2 — \/3/2, <PTv/i{ s ir/2) — ~V2/3, < / J "/ 2 ( s ^/2) 
•\/2 and a(s 7r /2,0) = \/2, which yields in this case 



n^+oo 2 5 / 4 \/7T 



[l + ^2A k (n/2)h k 



(2.22) 



fc=i 



Step 1.6. It remains to estimate the second term in (2.131, 



|A^ n + te je/3 )| 2 cft, 
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but a rough control of this remainder term will be enough. 

For n large enough, we can apply (2.7) uniformly with respect to t € [0,e] to get 

2 



Rn ' E ~ C|Mn| rnT 1 1 la +te'»/3|V2 



Ki„(|(-(^+*e*/»))»/» + ?[ 



C > , 



< C 

< c" 



e|Mn| 



1/2 



| _i + ee «/3H/2 



sm 



(| Mn |(-l + e e ^ 3 )) 3 / 2 + ^ 



C" > 



exp 



1 



ips(e) , C" > 



where we used 



sm -z 
3 



3/2 



< ^cxp(-?Rez 3 / 2 ) +o(l) 



for 27r/3 < |argz| < n. Thus, it exists C = C(0,e) > such that, for n large 
enough, 

R n ,e <Ch-^ 3 exp (-~<p e (e) 



Using that, for e < sg, 



(2.13) implies 



^ 3 /2 

-^fl(e) < -ipe(se) = g^V sin6>, 



\u n (t)\ 2 dt = I n>£ (l + (e- c /' 1 ")) 



for some c > 0. The expansion (2.21 ) finally yields the proposition: 
Proposition 2.5 Let u n (x) = Ai(fx n + xe l9 / 3 ) and h n — |/i„|~ 3 / 2 ; then 



1/6 



7 1— > + QO 



^exp — m 3/2 sin0 1 + ^^(0)^ 



k=l 



(2.23) 



This concludes the analysis of the numerator of (2.3). 
Step 2: The J n term 



It remains to estimate the denominator of (2.3). We first recover a real integral, 
using the analyticity of the Airy function and its decay at infinity in the sector 
| argz| < 7r/3: 



6/3 ul(x)dx 



+ 00 



e l9/3 Ai 2 {e ie/3 x + fi n )dx 



lim / Ai (z)dz, 
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where jr is the segment [fi n , fi n + e^^R] C C. The Airy function is holomorphic 
in C and 77?, can be decomposed into 

where Cr is the arc of a circle (centered at n„, of radius R) 

Cr : [0, 1] 3 t >-> M „ + Re lte '\ 



Hence, 



J8/3 



3t 



u 2 l (x)dx = lim 



+00 



Ai 2 (x)dx + I Ar(z)dz 



Ai (x)dx 



since, using (2.6), the integral J c Ai 2 (x)dx tends to as R — >• +00. 
Thus, 

f +00 



u n (x)dx 



Ai (x)dx, 



(2.24) 



since Ai(x) is real for real x. 

Now we use the formula of Lemma 2.3 which yields 

/•+00 

Ai 2 (x)dx = Ai' 2 (fi n ). 



(2.25) 



Here we have used that (2.6) and ( |2.8| imply 

xAi 2 {x)~ Ai' 2 {x) -> 

as x — >■ +00, and Ai(^„) = 0. 



According to the expansion (2.9) of Ai 1 , with the notation of (2.14), there exists 
a real sequence (d' k )k>o such that 



u 2 l (x)dx 



-/C 1/3 



3h n 4 



fe=0 



- sin 



3/i r , 



J EH^HiC 1 ■ (2-26) 



fc=0 



The expansion (2.10) of \fi n \ gives 



3tt 

n— >-+oo 2 



(n-1/4) [l + J2lk(n~l/i)- 2k , ( 7l ) t >ieM M (2.27) 



fc=i 



17 



Hence, for some sequence (Sk)k>i, 

(5c + 5)g- 1 )'<W?+-» 



cos 



o _\ +°° 



k=0 



(-1)" l + ^5 2fc (n-l/4)- 



2/v 



k=l 



From (2.26), we finally get 



Proposition 2.6 There exists (Bk)k>i such that 
u 2 n {x)dx 



n— >-+cxd 



lK^U + ^B k {n-llA)-A 



(2.28) 



where h n is defined by (2.1^). 



Step 3: Conclusion 

Writing each expansion in powers of h n a nd (n — 1/4) ~ 2 as expansions in powers 
of 7i -1 , we deduce from (2.3), (2.23) and (2.28) that for some (Cfc(#))fc>i, 



-\-oo 



-hi/ 2 



n exp 



3h n 



3/2 . , 

m a sin ( 



Finally, ( |2.27| ) yields 
4 



n eX P 



3/i„ 



3/2 . , 

m a sin ( 



— exp ^27rm^ 2 sin (n — 1/4) 
x(l + £c^)n- fe ), 



fe=i 



and the statement of Theorem |l.l| follows. The particular case of k^ ) (1,7t/2) 
follows from \2.22\ . 

The case of the Neumann realization is similar, putting h n — \v n \~ 3 / 2 and using 



(2. II) instead of (2.10). Notice however that (2.25) becomes in this case 

»+oo 

Ai 2 (x)dx = \v n \Ai 2 (v n ), 



hence we use (2.7) instead of (2.9) (this was wrongly stated in [IB]). 

In the next section, we will perform a similar analysis for the instability indices 
K n (2k, 9), k > 1. The asymptotic expansions of the Airy function will be replaced 
by WKB estimates on the eigenfunctions of A(2k,9). 
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3 Instability of even anharmonic oscillators 



We would like to understand the behavior as n — > +00 of the instability indices 



K n (2k, B),k> 1, using formule ( 1.6 ). In this purpose, we first consider the solutions 



Tph of the selfadjoint operator —h 2 4-^ + x 2k — 1. We will then explain the relation 



between tp^ and the eigenf unctions of A{2k,8) in Paragraph 3.2 



3.1 WKB estimates for the eigenfunctions 

We use the complex WKB method (see appendix) to understand the asymptotic 
behavior of the solutions iph € L 2 (M.) of equation 

V h (2k)1> h = (-h 2 -^+x 2k -l)ip h =0, (3.1) 

where k > 1. The eigenfunctions of the operator A(2k,9) are indeed related to 
the functions iph, since the transormat ion x H> | \^ 1 / 2k q 1 - 6 / ( 2k + 2 ) x maps formally 



A(2k, 9) into the operator Vh(2k) of (|3.l[), with h — |A| 2* (see the details in 



paragraph 3.2 1. Consequently, we are especially concerned with the behavior of 



ij)h{x) as x goes to infinity in the sector 

Se/(2k+2) = IxeC: |argx| < - }■ , (3.2) 



^/( 2 fe+2) - ^ t ^ • ^ 2 ^ fc - ^ 

but also at the limit h — > 0. 

Let = e^ n ' k , j = 0, . . . , (2/c — 1), denote the zeroes of the function a; H> 
x 2k — 1, and let us consider the function x n- \J x 2k — 1 defined on the set 

2fc-l 

£> = C\ (J e y,r /*[l,+cx>[. 
We then denote by S(z) the function 

/z 
y/ a; 2fc — 1 dx, 

where the integral is taken along a path defined on D joining 1 and z, the integral 
being independent of the path. 

If £ denotes the union of level sets (also referred as Stokes lines) of the function 
ReS 1 starting from the zeroes Xj = e lJ7r / fc , j = 0, . . . , (2k — 1), then C \ C consists 
of 2k + 2 unbounded connected components Sj , each domain X j being bounded 
at infinity by the asymptotic directions M+e l{2] - 1)v,{ - 2k+2) and IR+ e < ( 2 i+ 1 W( 2fc + 2 ) ) 
along with the two domains 

(2fc+i \ 
C\ |J E7j n{lmz^0} ; 
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see the following figures. 

Let r = D + UP_UE U E_i U E+i, and for e > 0, 



r e = {z G r : d(z, dT) > e}, (3.3) 
dT being the boundary of T. 

Then we have the following statement, where the asymptotics as \z\ — > +00 also 
hold for a potential x 2k — e v , (igl, instead of x 2k — 1. 



Proposition 3.1 



1. For all (ieR and rj G] - 7r/(2fc + 2),7r/(2fc + 2)[, i/^ sofoes pJp or 

rf 2 



(3.4) 



wi/i V'/i G -L 2 (R + e I? '), i/ien i/iere exists c G C suc/i i/iai, /or a// a G] — 
37r/(2fc + 2),37r/(2fc + 2)[, 

^(*) = ( Z M -1)1/4 ex P (~^ S (- Z )) + z = re ia , r^+c», (3.5) 

uniformly with respect to h. Moreover, ip^ is the unique solution satisfying 
(3.5) for a G] — it/ (2k + 2), -k / (2k + 2)[, and iph is holomorphic in C. 

2. In the case \x = 0, there exists a sequence of functions Uj defined on T, j > 1, 
such that 



iph(z) 



fs-s-0 (z 2k 
uniformly on T £ , for any e > 



^yyr/i fi + E^^j cx p (~^) 



(3.6) 
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Figure 5: Stokes lines of the anharmonic oscillator 



' dx 2 



X 



1. 



Proof: We apply Theorem A.l with f(x) = x 2k — 1, and with base points 
<Zo = 1 and ai = +ooe l,) the point at infinity in the direction arg -1 {?7}. The 
analysis of the Stokes lines show that for all z E S_i U So U Si, one can find 
a canonical path joining a\ to z : it is enough to paste a ray [i?e lJ) ,ai[, where 
R > is such that ReS(Re iri ) > ReS(z), with the path S^QS^), S(Re ir1 )]). 
This construction holds for the solutions of (3.1) as well as for those of (13. 4|, the 



asymptotic directions of the Stokes lines being the same in both cases. 
For /z = 0, we can also find a canonical path 7(2) joining a\ to any point z € 
2?_ U X>+, so that r C A (ao, a\) according to the notation of the appendix. Hence 
it remains to check that there exist k, p, M > such that, for any z £ V, the 
function a(z) defined in (A.3) satisfies (A. 4), and such that for all z e r e , there 



exists a canonical path 7(2) joining ai to z and satisfying (A. 5) 



Conditions (A.4) and (A.5) follow from the behavior of a as 



a(z) = Cz~ 2k - 2 (l +o(l)), CelR^O 



which concludes the proof of the proposition. 



□ 



The asymptotic behavior of the solutions of (3.4) being known in the whole 
sector arg -1 ] — ir /(2k + 2), ir/(2k + 2)[, we will be able in the next paragraph to 
justify the transformation mentioned above to switch from A(2k, 8) to r Ph(2k). 
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3.2 Another expression of the instability indices 

Let n > 1, A„ the n-th eigenvalue of A(2k, 9), and u n an associated eigenf unction. 
The scale change x = |A n | _1 / 2fc a; will first map the equation 

(A(2k,9)-X n )u n (x) = 

u n {h,- 



into 



" n dx 2 



-hi— + e l0 x 2k - e iargA 



0, 



where we put 



|An|" 



This will be our semi-classical parameter, observing that h n — > as n — > - 
Then we perform the analytic dilation y — e * 6 '/( 2fe + 2 )j j an d (3.7) becomes 

|A„|e^/( fe+1 ) {-h 2 n ^ + y 2k ~ e *(«^-fl/(*+i))^ = 0j 



(3.7) 

(3.8) 
-oo. 



where 



TP hn (y) = Un(h- 1/(k+1) e- i9 ^ k + 2 h)- (3.9) 

Sin ce u n € £ 2 (IR), we have ^h,, G L 2 (e _,,e ' ( 2fe+2 )R). According to Proposition 
we also have ^/j n G L 2 (R) , hence i/'/i„ is an cigcnfunction of the non-negative 



3.1 



selfadjoint operator 



K dy 2 



.2ft 



associated with the eigenvalue e l ( argA ™ e /( fe + 1 )). Therefore, we have necessarily 

e i(argA n -0/(fc+l)) = 1 



hence we recover the operator Vh{2k) defined in (3.1 1, and all the eigenvalues of 
A(2k,6) lie on the half-line arg^ 1 !^}. 

According to paragraph we can choose to normalize the eigenfunction u n so 
that the function ipf,, defined by (3.9) is characterized by (3.5) with c = 1, as 
x — > +oo for instance (or equivalently for \x\ — > oo along any ray included in the 
sector arg -1 ] — n/(2k + 2), n/ (2k + 2)[, see paragraph 3.1 ). 

The spectral projection associated with A„ b eing of rank 1 according to Lemma 5 of 

i 

[5J, Formula (1.6) holds for n n (2k, 9). Using (3.9) and the scale change 
we get 



K n (2k,0) 



k\^h n {e 2{k+1) v)\ 2 dy 



Now again we write the denominator as a real integral, as in Section 2.3 



g 2(k + l) 



^lAe^ TJ y)dy 



lim 



^hAz)dz 
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where 7^ is the segment [— i?e l2 < fc+1 ) , Re' 2 '- k + 1 > } C C. We write 

7fl - [-R,+R] vc B y(-c R ), 

where Cr(£) = Re lt ^+^ , t £ [0, 1], and the exponential decay of iph n as 
in the sector (3.2 1 yields 



-00 



lim 

R— ¥+00 



i> 2 h (z)dz = 



(3.10) 



(see Lemma 3.1 and notice that S(z) —> +00 as \z\ — > +00 in the whole sector 

((gi). 

Similarly, 



lim 

R->+oo 



Thus, we have 



Kn(2k,0) 



ipl {z)dz = 0. 



J R \*Ph n (e 2{k+1) x)\ 2 dx 



(3.11) 



Finally, the potential in A(2k, 9) is even, hence the solutions i/jh n are even or 
odd, and we can replace the integrals in (3.11 1 by integrals over K + . 
We summarize all the previous informations in the following proposition: 

Proposition 3.2 For all n > 1, the n-th eigenvalue X n of the operator A(2k,9) 
writes 

\ - r JS/(k+l) 

where r n > is the n-th eigenvalue of the selfadjoint anharmonic oscillator 

d 2 



,.2/v 



dx 2 



The following Weyl formula holds: there exists a real sequence (sj)j>i such that 

+00 



n— y +00 



l+J2*j(n + l/2)- 2j ) . (3.12) 
3=1 



There exists a unique eigenfunction u n of A(2k,0) associated with X n such that 
u n (x)=^ h M /(k+1) e m2k+2) x), 



where h n is defined by {3.8) and if) h is the unique solution of equation (3.1) satis- 
fying (3.5) with c = 1 as x — > +oo . Lastly, we have 



K n (2k,9) 



jg+ \^ h Ae W ^ k + 2 ) X )\ 2 dx 
Jr+ $h n (x) 2 dx 



(3.13) 



The expansion (3.12) is proven in |16) . Theorem (2 — 1), for operators of the form 
— ^27 + P( x )j where P is a real polynomial of even degree. 



In the next paragraph, we use the asymptotic expansions given in Lemma 3.1 



with h = h n — > 0, to get an asymptotic expansion for the numerator of (3.13). 
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3.3 Estimates on the norm of the eigenfunctions 

We assume without loss of generality that 9 > (if 9 < 0, replace 9 by \0\). 
For a fixed S > (which will be determined later i n th is paragraph), the hal f-lin e 
[S, +oo[e' i0 /( 2fe+2 ) belongs to the set T e defined in < 3.3 ) , e < S, hence using (3.6|, 
we have 



+ exp(-is(x) 



(3.14) 



uniformly for x € [<5, +oo[e j6, /( 2fe+2 '. 
Then we write 



+oo 



\il>h(e l *vxVx)\ 2 dx = h{h) + R s (h) 



(3.15) 



where 



h{h) 



+ 00 



\^ h (e^("+^x)\ 2 dx, Rg(h) = / |V>h(e i2 < fc+1 > ir)| 2 (fa: 



and we first estimate Ig(h). The expansion being uniform with respect to x, we 
can take the integral over [S, +oo[ in (3.14). Thus there exists a sequence (vj)j>i 
of functions such that 



h(h) = / ag(x,h)e^ ipe - k{x) dx 



(3.16) 



where 



and 



a e (x, h) 



f 1 + V vj (x)h j 

& "1| 1/2 I 



(fk _ 1 j 1 /2 tft< 



Now we apply the Laplace method to Is(h), checking first that the phase (fg^ has 
a unique non-degenerate critical point in [5, +oo [. 
We have 

WA*) = -£ - !| 1/2 cos Q arg( S 2fe e^ - 1) + yj^) 



hence 



Ve,k{x) = -l^e 1 ^ - II 1 / 2 cos Q arg(x 2 V & - 1) + ^^yy) 
since arg , 7r] , we have 



- arg(x e — lH : r 

2 8V 7 2(fc + l) 



6> 7T 



2'2 2(fc + l) 
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and ip'g k (x) vanishes if and only if 



-arg(. e W -l) + _ 



7T 

2' 



that is 



x 2k sm(k9/{k + 1)) 
x 2k cos(fc0/(fc + 1)) - 1 



= — tan ■ 



which gives a unique critical point x@ y k for ipe.k, 

tan(0/(fc + 1)) 



sin(fc0/(fc + 1)) + cos(k6/(k + 1)) tan(0/(fc + 1)) 



(3.17) 



(3.18) 



One can easily check that this is a non-degenerate maximum. 

Of course ipg^(5) < ipe,k( x 6,k) if 6 < xg.k, an d for M > xg.k, the integral 

Jm°° a s (x, h)ei ve ' k{ - x > h) dx can be est imate d as in (2.20). 

Thus, the Laplace method applies to (3.16), and there exists a sequence (rj(2k, 0))j>i 
such that 



h{h) 



2?r 



-f- OG 



h^o \(xl k k e^/(k + i)_ 1)ip >> 



(3.19) 



The following lemma gives a rough estimate on the remainder term R$(h) in 



(3.15), provided that 5 is chosen small enough. 
Lemma 3.3 There exist S > and c s]0, 2tpg,k{xe,k)[ such that 

Rs(h) = 0{e c ' h ). 



(3.20) 



Proof: The function ip^ being holomorphic in C (see Proposition 3.1), we can 
apply the maximum principle in a complex neighborhood of [—1, 1]. 
Let a > 0, and 6 > small enough so that 



Vx € [-1 - a, 1 + a] + i[-5,5], ReS(x) > ReS(xg^e 2 ^ 



(3.21) 



Let fl C C be the open ellipse centered at 0, with axes [—1 — a, 1 + a] and [—i6, iS]. 
Let "fd = dfl D T e , where T e is the domain defined in (3.3), and j g = dQ \ jd- 
According to Proposition 3.1 we have 

1 



+ 00 



Va; G 7 d , tfj h (x) 



h^O (x 2k - 1)V4 



1 + / ]uj (x )hP exp 



-S(x) 



Hence, using (3.21 ), 



Va; € jd, \ip h (x)\ 2 = 0(e Cd/h ) where c d < 2y e>k {xe,k)- 
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On the other hand, the symmetry of cigenfunctions yields 

Vx£j g , \iph(x)\ 2 = 0{e c » /h ) where c g < 2<pe :k (x e>k ), 
and the maximum principle applied to iph in Q gives 

sup \Mx)\ 2 < sup \Mx)\ 2 = 0{e c ' h ), c < 2tpe, k (x e ,k), 

x£Q xGdCl 



which yields (3.201 as soon as 5 is small enough to ensure that Z?(0, S) C Q. B 




Figure 6: The open set SI. 



Thus (3.151 and (3.19) lead to the 



Proposition 3.4 There exists a sequence (ry(2fe, #))j>i such that 
/ IMe^^x^dx ~ C fc (0)/i 1/2 e^ \1 + Yr ] (2k,e)h j 



where 



and 



Ck(9) 2 \{xf k ^/i^)-l)^ k {x e<k )\y 



(3.22) 



c k {6) = 2(pg ik (xg ik ). 

In the following paragraph, we get an asymptotic expansion for the denomina- 
tor of (l3~13|). 
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3.4 Solutions of the selfadjoint oscillator on the real axis 



The WKB expansions established in Proposition 3.1 are no longer available on 
the real axis, because (3.6) does not hold in the neighborhood of th e Stokes line 
[— 1, 1]. In order to estimate the denominator L °° tp^(x)dx of (3.13), we will use 
another representation of the function iph, like in [T3] (Exercise 12.3). 

Lemma 3.5 There exists 8 > such that, as n — > +oo, 



(3.23) 



where the function u(-, h) satisfies: 



• u(-,h n ) is an asymptotic solution of (3.1), namely 

V hn (2k)u{-,h n )=0(h™). 

• There exist C°° functions ip± defined in a neighborhood of such that 

1 



(%/>±(0) 2fc + £ 2 -i = o, ^±(o) = o, 

and functions 



5^-(o) 



-1 



(3.24) 



+ 00 



such that 

Vx e [±1-25, ±1 + 25], u{x,h) 



'2%h 



e i(^+Mt))b ± {£,h)d£. (3.25) 



• There exists C > and a sequence of real functions (aj(x))j>o, where 

a {x) 



(3.26) 



such that, for all x G [— 1 + 5, 1 — S], 

/ 1 \ + °° 

u(x,h) = cos[-<p(x)-^)J2(-l) P a2p(x)h 2p 

\ ' p=0 

)+oo 



•sin [ ^¥>W 



where 



/x 
VI - £ 2fc di. 
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• For all x < —1 — 8 and all x > 1 + 8 , we have 

u(x,h) = 0(h°°). 

Proof: The construction of u(-,h) is suggested in Exercise 12.3 of [13] in the 
general case of a simple-well potential. By construction, 

V h (2k)u(;h) = 0(h°°), 

and the spectral theorem implies that there exists a(h) such that 

\\ij h -a(h)u(;h)\\ LHm =0(h°°). 



(3.28) 



To show (3.231, it is then enough to check that a(h) = 1, which can be made 
by connecting the two solutions tph and u(-,h), comparing them both to a third 
solution expressed in terms of Airy function (see [21], Chapter 13, paragraph 7). 
Let us recall that iph is the solution of Vh(2k)tph = characterized by its decay at 
infinity in the sector arg -1 ] — it/ (2k + 2),7r/(2fe + 2)[. On the other hand, u(-, h) 
is an asymptotic solution such that, for x G [— 1 + 8, 1 — 8], 



u ^ h ) = (1 _ %y/i cos 



\-\ip{x))-Oih). 



(3.29) 



Another solution Wh is given by the following expression in the neighborhood, say, 
of the zero +1 (see [2"T]): 



(x) = ( T 



1/4 



Ai 



where 



-,S(x, 



(cm 
U 2/3 

2/3 



0(K 



(3.30) 



Expression (3.30) actually holds for all x G [— 1 + 8, +oo[, and comparing it to 
(|3.29|) as h 0, we get 



to/, (as) 



(1 + 0{h))u{x, h). 



On the other hand, comparing (3.301 with (3.5) as x — > +oo, 

w h {x) = ^{l + 0{h))^ h {x). 
Thus, comparing the principal parts, 



iph(x) 



C 



u(x,h)+0(h°°), 



that is the coefficient a(h) in (3.28) is independent of h. Up to renormalization of 
the function b(x,h) defining u(-,h), we can assume that a(h) — 1 in (3.26), and 
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(3.23) follows. 



This lemma enables us to compute the norm of tph on the real axis, using the 
expressions of u(-,h). Since expressions (3.25) and (3.27) differ from a 0(h°°) 
term on the intersection of their domain of validity, we can write 

u = u_x- + UoXo + u+x+ + 0(h°°), 

where (x-tXo,X+) is partition of unity and u±, uq are respectively the functions 
on the right-hand sides of (3.25) and (3.27). Thus, 

\\u(-,h)\\ 2 = I \u_ X -\ 2 dx+ I \u Xo\ 2 dx+ f \u +X +\ 2 dx 



+2 / \u \' X 0X-dx + 2 / \u^x x+dx + O{h™) 

I-(h) + I (h) + I+(h) + I 0t _(h) + / ,+W + OQT). (3.31) 



In order to compute J_(/i), we use (3.25) to write the norm as a triple integral as 
in 0: 



I-(h) 



1 

2nh 
1 

2~Trh 



w_ (x)ii— (x)x- (x)dx 



e i(^(«-e')+V'(?)-^(e')) 



[-l T 5] x xM ( e,0 



x&(f , h)b(£, h)x-(x)dxd£d£, 
hence after the change of variables (dependent on £) 

i = e-e 

X = 35 + 9^(0) 



I-(h) = 



1 

2nh 



e^ £ &A(x,i,£,h)dxdn d£, (3.32) 



where 
and 



A(x, I £, h) = b% h)b(i + h)xl(Z - %K0)- 



We now apply the stationary phase method to the integral with respect to (x, £). 
First we have 

ds^ (x, = -| and f) = -5 + %/>(£) - 9^(1 + 0. 

hence has a unique critical point at £) = (0,0). 
Moreover, for all (ft. 



Hess $^(0,0) = 







-l -aim 
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hence (0, 0) is a non-degenerate critical point, and the stationary phase method 
applies, with $^(0,0) = 0. 

Thus, it exists a sequence (mj(£))j>i of functions such that 

J_(fc) ~ / \b (0\ 2 \x-(-d^(0)\ 2 ll + ^mjWhAdt, 

^° *'» y j=i / 

that is 

I_{h) ~ ^cL^, c° t^O. (3.33) 

Similarly, we show that 

J +W,~ E4^. 4^0. (3.34) 
j=o 



We now handle the integral Io(h). According to (3.27), there exist real se- 
quences (a'j(x))j>o, (a'!(x))j>i and (o-j{x))j>i, where 

C 2 

a' {x) 



VT— x 2/c ' 

such that 



Uo(:r,fc)| a ft ~ Q cos 2 f-^i) - J) 



2 > v 



cos ( -^(i) ) E 5 j( a; )' l2j+1 - ( 3 - 35 ) 



Let us integrate the first term of the right-hand side, which contains the principal 
part of the expansion of Io(h). 

I xl{x)cos 2 (hp(x) - ^) Y^a'^h^dx 

»~ \ L xl(x) (' + * l v[x) ) IX-m* 2 ^ 



3 
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where 



If, 2 , +00 

R(h) = - / Xa( x ) a '( x > h) sm T<P( x )dx, a'(x,h) ~ > a'j(x)h J 
2 J R h h ^ j~ 

Successive integrations by parts (principle of non-stationary phase) yield 

R(h) = 0(h°°), 

hence 



f /l 7r\ +CX) 

y R Xo(z) cos2 (^te) _ I) a '( x > /l ) da; h ~ n i 



h 2 3 



for some real sequence {fij)j>o- 

Performing the same computations on the two other terms in (3.35 1 finally yields, 
for some sequence (c? )j>o, 



Lastly, this method gives the same kind of asymptotic expansions for integrals 

C 2 Xo{x)x±{x) 



(3.36) 



Jo,±(fe) in (plj): 



+00 . 



2k 



-dx. 



(3.37) 



Gathering (|3.33|), (|3.34|), (|3.36|) and (|3.37|), according to (3.31 ), there exists (cj)j>o, 

(3.38) 



c 7^ 0, such that 



+00 



\H;h)\\ 2 ^Y,^ 



Thus, using (|3.23|), we have proved that 

p+00 



/+00 +0° 
\ip h (x)\ 2 dx ~ ^cjh 3 , c 7^ 



(3.39) 



3.5 Conclusion: Proof of Theorem 11.21 



Expressions (3.11), (3.22) and (3.39) yield the following statement: 

Theorem 3.6 Let k e N* and 9 satisfying < |0| < (k + l)ir/2k. If n n (2k, 
denotes the norm of th n-th spectral projection of operator 



A(2k,9) = -^+e* 9 x 2k 
dx £ 
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and X n its n-th eigenvalue, then there exists a sequence (C'-(2fc, 0))j>o such that 



+00 



K„(2k,e) ~ hl l /2 e c '^ /h "YCJ2k,9)hi, C' Q (2k,6)>0, (3.40) 

3=0 



where h n — \X n \ ^ and 



where 



xe.k 



tan(0/(fc + 1)) 
sin(fc6>/(fc + 1)) + cos(fc6»/(fc + 1)) tan(0/(A; + 1)) 



£ e 2(fc+l) 

Im / (l-* 2fc ) 1/2 <&- 
Jo 



(3.41) 

(3.42) 
(3.43) 



Finally, we get Theorem 3.13 by replacing the asymptotic expansion (3.40) by 
an expansion in powers of n~ L , using the following quantization rule: 

Lemma 3.7 For all k > I, there exists a sequence (s J k )j>i such that 

+00 



1 

h„ n—t+oo 



{k + 1)V ^^(n + l/2))[l + Y t 4(n + l/2) 

3=1 



-23 



This formula follows from (3.8) and (3.12) 



In the last section, we prove Theorem |1.4| and Corollary |1.5 



4 Completeness and semigroups 
4.1 Completeness of eigenfuctions 



In this paragraph we prove Theorem |1.4| First of all, let us recall that, if H is 
an Hilbert space and p > 1, the Schatten class C p {%) denotes the set of compact 
operators A such that 



'+00 



l/p 



Mil 



J2lJ-n(A) P ) <+OC, 



(4.1) 



where ((i n (A)) n >i are the eigenvalues of (A*A) 1 ^ 2 , repeated according to their 
multiplicity (see [E]). The space C P (H), p > 1, is a Banach space. One can also 
define the space C P (H) for p e]0, 1[, but the quantity (4.1 1 does not define a norm 
anymore, and C P (T-L) is not a Banach space. 

We already know that the resolvent A(2k,8)^ 1 is compact for any k > 1 and 8 
satisying (1.5). We now prove like in [52] that it actually belongs to a Schatten 
class: 
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Lemma 4.1 For any e > 0, \9\ < and fc > 1, we ftaue 

(X^M)) -1 e C^ +e (L 2 (R)). 

fc+i | - 

Proof: Let us show that, for all e > 0, the series /J,„ 2t: is convergent, where 
(A*n)n>i ar e the eigenvalues of 

([(^fc,^))- 1 ]*^^))- 1 ) 172 = ([^^(^fc,^))*]- 1 ) 172 . 

If {v n )n>i denote the eigenvalues of A(2k, 6)(A(2k, 9))*, then we have to check 
that 

£ -« p/2 < +°° 

n=l 

as soon as p > . 

.A(2fe, 6)(A(2k, 6))* is a selfadjoint operator, and if p(x, £) denotes its symbol, we 
define its quasi-homogeneous principal symbol -P(x, £) as 

following [23 . 
Then we have 

P(x,0 = \e + e w x 2k \ 2 = e + 2cos6ex 2k + x ik , 

P{r 1/4k x, r 1/4 = rP(x, f), T > 0. (4.2) 

Moreover P is globally elliptic, in the sense that 

V(x,O^(0,0), |P(a;,OI>0. (4.3) 

Hence the results of [53] allow us to apply the following Weyl formula: 

N{t) := #{j > 1 : Vj < t} ~ / 



which, with £ = j/„ and using (4.2), yields 

fe+i 

n ~ Cv n ik 
n— >-\-oo 

where C= Vol P-^O,!]). 

Thus the series ^2 v n P ^ 2 converges if and only if 

2kp 



n k + 1 < +oo, 



n=l 



that is if and only if p > . 
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Since the operator A(2k, 9) is sectorial and its numerical range is included in 
the sector Sg = arg _1 [0,#], the resolvent estimate 

\\(A(2k,B)-X)- 1 \\=0(\X\- 1 ) (4.4) 

holds outside Sg, and if we denote p — + e, then 

(k + l)n it 7r 

6 < 2k < ^T7 = v (4 ' 5) 

for e small enough, as soon as k > 1. 

Consequently, Theorem 1.4 follows from Lemma 4.1 and Corollary 31 of [TU], p. 
1115. 



In the next paragraph, we prove Corollary |1.5| 

4.2 Semigroup decomposition and control of the remainder 

Theorems 1 1 . 1 1 and 1 1 . 2| for m — 2k, with k = 1/2 or k > 1, yield as n — s- +oo 



- tA "n„|| = e- mcX - Kn (m,6) 



= ^K^) e _ tR eA„ e c fe ( 9 )|A n |* (1 + o(1)) 

= ^L^e-^l^^e^WIAnl^^ + o^)) ( 4. 6) 

where, in the case m = 1, the constants 0) = and Ci/2(0) = C(0) are 

those of Theorem 11. II 

Thus, for k = | and for any f > 0, there exists _/V t £ N such that 

Vn>N t , e -*|A»|cos^ e c 1/2 (0)|A„|§ > e ^|A„|^ 



and (4.6) is the general term of a strongly divergent series. 



For k = 1, yields 

||e- tA »n„|| = ^Ml e ( c 1 W^ c ° s f)l A "l(l + (l)), 

hence the series ^(i) is norm- convergent t > T and not norm-convergent for 
t < T. 

Lastly, if k > 2, there exists N t such that 



and (4.6) is controled by the general term of a convergent series. 
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To check that the series £ m (t) (when convergent) converges towards the semi- 
group associated with A{m,6), we use the density of the family (u n ), where the 
cigenfunctions u n are assumed to be normalized by the condition (u n , u n ) = 1, so 
that (u n ,u n ) n >i is a biorthogonal family (see [6]), namely 

Vn, m G N, (u„, u m ) = 6 n>m . (4.7) 

Then we have 

e -tA(mM) Un = e -tX nUn 

and on the other hand, 

E TO (t)ii n = y^e^^IIj-Un = e~ tXn u n . 
j=i 

Here we used the formula 

(see [B], [1]) which holds for rank 1 spectral projections, together with the biorthog- 
onal property (4.7). 

Hence by linearity, e - tA ( m -fi) and £ m (i) coincide on Vect{u n : rt > 1}, and hence 
on T>(A(m,9)) by density (see theorem |l.4[ ). 

To prove the inequality on the remainder, we write e~ tA ( m ' S '(I — H<n) as the 
remainder of order N of the series E m (t), which we estimate roughly to get (1.151. 



A The complex WKB method 

Here we recall the main theorem used in section [3] The general theory of complex 
WKB method is detailed in [3T] , [23] , [T2] , and we will follow the point of view of 

Let us first recall the following convention: for (j G K and x G C, a path joining 
+ooe l13 to x will denote a path 7 :] — 00, to] —> C of the complex plane such that 
j(to) = x and such that 7Q — 00, T]) coincide with the ray [R, +oo[e z/3 , for some 
T <t and some R > 0. 

We study the asymptotic behavior as |x| — > +00 or h — > of the solutions iph of 
equation 

-h 2 ^n{x) + f{x)^ h {x) = ^ xeC, (A.l) 

where / is an holomorphic function on C. Let T> be a simply connected subdomain 
of C \ / _1 (0), and consider a determination of the function y/J on I?. For a fixed 
ao e P (or ao a zero of /), we put 

S(z) := T V?(t) (A.2) 

where the integral is taken along any path defined on T>. 

Let ai € C, with possibly a% = +ooe l/3 , j3 G K. We call canonical path joining ai 



3G 



to z, and we denote by 7(01, z), a path of the complex plane joining a\ to z such 
that the function ReSo 7(01, z) is decreasing. We then denote 



and for e > 0, 



A(ao, ai) = {z G C : 3 7(01, z) canonical path }, 



A e (a ,ai) = {z e A(o ,Oi) : <9A(a , ax)) > e}, 



where 9A(ao,ai) denotes the boundary of A(oo,Oi). 
We also define the function 



1 



/3 /4 



1 



fl/A 



(A.3) 



Then the following theorem holds, see [3T], Theorem 11.1 (p. 222) and Chapter 
10. 

Theorem A.l (WKB estimate) Let f, ao, a\, A(ao,ai) and A £ (ao, a\) be as 
above, and S C A e (ao,<n) (not necessarily bounded). If ax = +ooe l/9 , we assume 
that Re S(ye zl3 ) — > +00 as y — > +00. We also assume that there exist positive 
constants k, p and M such that, for all z G A(a Q ,ai), 



k 



and for all x G S, 



i + |s(z)| 1+ '' 

S'(x)\ 



+p 



dx < M, 



(A.4) 
(A.5) 



where 7(01, z) is a canonical path joining a\ to z. 

Then, there exists an analytic solution tph in C of equation (A.l) such that: 
1. For all a G K satisfying 

3R a >0, [R a ,+oo[e la C E, 

we /lave 



1 



exp 



1 



S(*) (l + o(l)), 



(A.6) 



uniformly with respect to h, as z goes to infinity in direction arg _1 {a}. 
Mo reove r, ifReS(y e ) —> +00 as y — > +00, then iph is the unique solution 
of (A.l) satisfying (A.6), and for any solution tph € L 2 (M. + e ta ) ; there exists 
c G C such that iph = ciph- 

2. There exist functions Uj = Uj(ao,ai), j > 1, such that 

+00 



1 



uniformly with respect to z in E 



1 + y^ j Uj{z)W I exp 
3=1 



-S(z) 



(A.7) 
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